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Abstract 



Conservation of baryon number as topological charge in effective chiral field 
theories imposes a local constraint on the time evolution of field configurations 
in the commonly used 0(4) model. Possible consequences for the formation 
of chiral condensates after a quench from random initial configurations are 
discussed. It is argued that efficient dissipative terms are necessary but not 
sufficient to unwind the randomly curled-up initial configuration, before col- 
lective motion towards the condensate can grow. The existence of soliton 
CN \ stabilizing mechanisms will further prolong or prevent this process. 



°V 

An attractive feature of the nonlinear sigma model and its extensions to higher chiral 
orders [[j] is the possibility to identify baryon number with the winding number [[| that 
Q_i! characterizes the map of compactified coordinate space onto the £77(2) manifold of the 
chiral fields $ = (a, 7r), as long as they are constrained to the 3-sphere a 2 + it 2 = f 2 . This 
concept has a profound basis |§ in the anomaly structure of the underlying QCD, and it 
has led to numerous successful applications as efficient description for meson-baryon systems 
without explicit fermionic fields. 

In the 0(4) linear a-model H, on the other hand, the fields (a, 7r) can freely explore the 
full 4-dimensional chiral space which allows for a convenient description of the possibility to 
restore the spontaneously broken chiral symmetry at finite temperatures to study critical 
fluctuations of the chiral condensate near the phase transition, or to follow the dynamical 
formation of the condensate after a quench. 

Unfortunately, with the unconstrained embedding of the chiral field $ into an R A man- 
ifold the concept of baryon number B as topological index is lost because the topological 
connectedness of this manifold is trivial. Of course, in a parametrization like 

$ = R(x, t)U(x, t), U e £77(2) (1) 

one still can maintain the definition 
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B= 2^I eijktrL i L i L * L > = ( 2 ) 

but this 'baryon number' B changes its value whenever the field R(x, t) as a function of time 
passes through the origin R=0 for some value of x. 

This is easily visualized in the 0(2) model in 1 + 1 space-time dimensions with two chiral 
fields, $ = (a, 7f), 

a(x, t) — r(x, t) cos cj)(x,t), 7r(x, t) — r(x, t) sxn<p(x,t). (3) 

All field configurations with boundary conditions 

a(—oo,t)—a(+oo,t) — f 7r , 7r(— oo, t) = 7r(+oo, t) = 0, (4) 

represent closed loops embedded in the a-ir plane which can be contracted continuously into 
the vacuum configuration a = f n , it = 0, and the 'baryon number' 

1 r°° d 

B = * Le&*- t)dx (5) 

changes by one unit whenever the sling is pulled across the origin. 

Soliton and vacuum configuration may be separated by a potential barrier which could 
prevent for classical dynamics such baryon number violating processes. Cohen |J has argued 
that the corresponding quantum mechanical tunneling amplitude vanishes in tree approx- 
imation in the limit where the number of colors N c —>■ oo. But a possible suppression for 
finite N c through quantum corrections to our knowledge has never been proven. By suitable 
changes in the parameters of the model (e.g. by increasing the temperature) the process can 
become even classically allowed. This can be nicely demonstrated in the 1+1- dimensional 
0(2) model with a temperature dependent value of f%(T) because in this model the stabi- 
lization mechanism for solitons breaks down already at temperatures which are much lower 
than the critical value for restoration of chiral symmetry. The lagrangian is 

£($) = ■ - V($) (6) 

with the potential 

V{^) = -^ 2 -ff-Ha. (7) 
In order to keep the minimum of V for finite H at <3> 2 = f%(T) we define 

In 1+1 dimensions this model stabilizes solitons, as long as the inequality 

jji < 0-047 (9) 
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is satisfied. They correspond to bound trajectories of a classical point particle moving in 
the potential —V, starting from the maximum of —V and returning to it after an infinite 
time. These configurations sling around the hat of the potential V close to (but inside of) 
the bottom of the valley. In terms of a- and 7r-masses 

< = T> ml = 2\fl + ml (10) 

we find from (|9]) approximately m a > 6.5 m n , as stability condition. The angular winding 
length L w ~ (<f)'(x = 0)) _1 (i.e. the inverse of the gradient of the angular field near the 
center of the soliton) defines the typical spatial extent of the stable 'baryon'. For fixed ratio 
m a jm^ it scales like (f^/H) 1 ^ 2 . In the symmetry limit H — > 4>{x) approaches tt(x/L + 1), 
so in this limit the soliton completely occupies the spatial box of length 2L. With increasing 
L w the radial bag r(x) gets very shallow, i.e. r(x) stays very close to f n (cf. fig.l). 
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FIG. 1. Soliton profiles in the 0(2) model (m CT = 1000, m n = 138, fo = 93) for three different 
values of /„- = q/o ( full lines: a=1.0; dashed lines: a=3.0; dotted lines: q=0.5 ). In the lower half 
the radial profiles r(x)/f, in the upper half the corresponding angular profiles 4>(x)/ir are shown. 

Suppose that H, A and /„- satisfy condition (H) for T = and suppose we have a 'baryon' 
represented by a stable trajectory with B = 1. With increasing temperature T, the inclusion 
of loop corrections into a renormalized C can lead to a decrease in f%(T) while the coupling 
constant A and H remain unrenormalized. This causes the soliton radius to decrease. At a 
certain temperature the value of H/(\f%) will exceed the critical stabilizing value @, the 
sling collapses to the vacuum point a = f n ,n = 0, and the 'baryon' has disappeared. 
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In the 3+1-dimensional 0(4)-model soliton stabilization requires terms of higher chiral 
order, like the fourth-order Skyrme term || 

C^ = ^jtr[L„L u rd 3 x. (11) 

Here the spatial extent of the angular chiral soliton profile is determined by (f n (T)e)~ l 
where the Skyrme parameter e remains essentially unchanged with increasing temperature. 
Thus the Skyrme stabilization mechanism works in the opposite way as compared to the 
stabilization through the pion mass in the 0(2) model discussed above: The angular winding 
length L W (T) ~ f^(T)~ l as the characteristic length scale for stable field configurations now 
increases with increasing temperature. But as before, the radial bag with increasing L w gets 
very shallow, i.e. R(x) stays very close to f^(T) as f n (T) goes to zero. So for the static soliton 
the violation of baryon number occurs together with the restoration of chiral symmetry, but 
classical fluctuations R(x, t) around the static configuration will cause fluctuations in baryon 
number already at lower temperatures. 

Of course we do not want baryon number to be violated just by increasing the temper- 
ature. One possibility to avoid this would be to strictly stay within the nonlinear sigma 
model, i.e. to impose the constraint a 2 + 7r 2 = f 2 (T) for all configurations at all temper- 
atures J7|. It is, however, not clear how to describe in such a framework the system in the 
phase where chiral symmetry is restored. 

A natural way to combine the convenient features of the linear a-model with the topo- 
logical advantage of the nonlinear a- model is to insist also in the 0(4) model on the angular 
and polar nature of the chiral fields, i.e. on the parametrization ([!]). Because for R=0 the 
angular field U is not defined the field configuration $=0 then is excluded from the man- 
ifold of possible chiral fields. This restores the nontrivial topological connectedness of the 
nonlinear cr-model. In our simple example of the 1+1 dimensional 0(2) model this means 
that we obtain topologically distinct classes of field configurations and consider the discrete 
manifold 

r = f r , = n2vr, (12) 

(n integer) as a set of physically distinct degenerate vacua. The proper space in which the 
chiral fields then live can be visualized either in polar representation (r, 0) as the half-plane 
(r>0, — oo < 4> < oo), or, in cartesian representation (a, n) as multiple sheets tied 
together in a multiple branch point at the origin. 

If we now repeat the above considerations, starting with a stable baryon at T = 0, 
increase the temperature, then at the critical value (|9]) for f n (T) the sling now collapses to 
a narrow hairpin around the origin which connects the vacua (|i~2"D with n = and n = 1, 

(j)(x) = 2tt0(x); 

r(|x| 0) = +0, r(x^±oo) = f n - Ae'" 1 "^ . (13) 

The stable nontrivial solution for r(x) represents the 'bag' profile of the remaining baryon 
with B = 1. Although its angular winding length L w has shrunk to zero, the radial bag 
still has a nonzero radius in coordinate space which increases with temperature as long as 
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the cr-mass decreases (in this particular stabilization mechanism which is peculiar to 1 space 
dimension). In fig.l soliton profiles are shown for some values of f%(T) above and below the 
critical value. The stable soliton solutions above the critical temperature do not exist if the 
point cr=7r=0 is not excluded. 

Of course, the exclusion of R=0 has severe consequences for the dynamical evolution 
of classical configurations. Recently it has been suggested that the exponential growth of 
collective amplitudes in the dynamical evolution of chiral field configurations after a quench 



may lead to disoriented chiral condensates in macroscopic regions of space [T5|]. It appears 
that these calculations have been performed [|rj]] and confirmed |T5| in the framework of the 
linear sigma 0(4) model Q i.e. with chiral fields embedded in a simple R 4 internal space 
without any constraints. It is an interesting question to ask whether and how the topological 
aspect discussed above may affect the results of such calculations. 

The typical difficulty appears already in the formulation of initial conditions: Commonly 
they are taken as random Gaussian ensembles centered around the origin <j=7r=0 to reflect 
the chiral symmetry restored at temperatures higher than the chiral phase transition tem- 



perature ||14| . In an (r, <p) representation of the chiral fields this would correspond to a 
uniform deviate in the angle (p. It remains, however, ambiguous over how many multiples 
of 27r this uniform deviate extends. Chiral symmetry would require to allow for all angles 
— oo < <fi < +oo with equal probability, which would imply that the initial field configuration 
could curl arbitrarily often back and forth around the origin from one lattice point to the 
next, with only the net winding number n of the configuration inside the box of size 2L, i.e. 
the difference 27m = <p(L) — 0(— L), fixed. Evidently, this is not a very meaningful choice 
for initial configurations. 

The high value of the temperature T reflected in the initial configuration implies a small 
correlation length £ ~ T^ 1 which characterizes the correlation function at t = 

< e i(p(xo+x,0) e -i<j}(x ,0) >= e -\<[4>(x +xfl)-<t>{x fl)] 2 > _ e -|a;|/£o_ 

This shows that the mean square net winding increases linearly with x. On a lattice with 
lattice constant a this corresponds to a random walk with \x\/a steps with a mean square 
step size of 2a/£o- Identifying £ with the lattice constant we therefore sample initial con- 
figurations for the angles (f>(x, 0) by random walks with Gaussian deviate for mean square 
stepsizes of order two, and fixed value of total net winding n. If a << 2L these config- 
urations again wind many times back and forth around the origin, but most of the phase 
differences between neighbouring lattice points are less than 27T, i.e. in other words, the 
local baryon density (~ 4>'(x, 0)/(2ir)) is small. 

Starting from such initial configurations, the singular point at the origin causes a severe 
local barrier against the free evolution of the classical trajectories: the growth of collective 
motion in radial direction towards the condensate can only take place (even for n = 0) after 
unwinding the random multiple twists around the origin. As long as the field is curled up 
around r = 0, high- and low-/c 2 Fourier components will have comparable power. This can 
be verified in the 1+1 dimensional model discussed above, if the time evolution is performed 
on a lattice in (r, 0) representation. The actual value of the total winding number n of the 
trajectory is of minor importance for the general features of the evolution during this stage: 
as long as the field curls many times back and forth around the origin, it is not essential at 
which multiple of 2tt it ends. 
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FIG. 2. The time evolution of the angular field in the 0(2) model from a random initial con- 
figuration (p(x, t = 0) (dotted line) to a quasistable soliton-antisoliton configuration 4>(x, t = 100) 
(full line) for strong symmetry breaking (m a = 2000,7717,- = 1000, 7^ = 93). 



The models discussed above do not even contain a term which would drive the unwinding 
of the curled-up configuration. The evolution of the radial part r(x) towards the condensate 
is driven by the slope of the potential (0) whereas there is no corresponding driving force 
to unwind the random large amplitude fluctuations in the angle 0(x, 0). This is easily 
demonstrated if (for the moment) we keep the radial field fixed at some value r(x) = r and 
then separately consider the equation of motion for the angular part 

= 0"-— sin (15) 

which is of course the sine-Gordon equation and, as a typical wave equation conserves the 
mean amplitudes of randomly fluctuating angles. 

Even the addition of a dissipative term on the left hand side of ( |i~5"l) (which could be 
linked to pion emission processes) to enforce the annihilation of large narrow fluctuations 
does not necessarily lead to the unwinding of the trajectories if the underlying lagrangian 
stabilizes solitons within the box of length 2L. The dissipative term will relax the random 
fluctuations efficiently into a quasistable ensemble of n + solitons and n_ antisolitons (with 
n = n + — nJ) distributed randomly within the spatial box. Although the configuration then 
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is characterized by a phase correlation length of order L w it generally will contain many 
sections winding back and forth around the origin which still prevent a definite orientation 
of the chiral field (even for n — 0). Whether optimal unwinding (n + = n) can be achieved 
depends on the ratio of three different length scales: the initial phase correlation length a 
(~ the lattice constant), the winding length L w of the stable (anti)solitons, and the size 2L 
of the box. A necessary condition for unwinding is a << L w , because otherwise the initial 
random walk configuration will be close to some quasistable soliton-antisoliton trajectory 
with a comparable large mean square deviation L/a. If a « L w « L the quasistable 
trajectory can be very different from the initial incoherent fluctuations but still may contain 
many multiply curled-up sections. Further unwinding due to soliton-antisoliton annihilation 
will only take place if at some point a soliton has sufficient overlap with an antisoliton. Only 
if L w is comparable or larger than 2L/ (n+2) then no stable soliton-antisoliton configurations 
fit into the box and the dissipative term can completely unwind all fluctuations. As we have 
mentioned earlier this is always the case in the symmetry limit (H — 0). 



1200 



1000 



800 



600 



400 



200 




FIG. 3. The power of some of the 256 Fourier components \itk\ 2 (k=2, 11,21,31) of the pion 
field 7r(x, t) (3) (from t = up to t = 15) for the same evolution as shown in flg.2. 

Fig. 2 shows a typical initial configuration for the angle <fi(x, 0) with n = for 256 lattice 
points. At the same time r(x, 0) > is chosen as narrow Gaussian deviate around r = 0. 
Without dissipative terms 4>(x, t) at much later times still looks like another random walk 
with similar characteristics, i.e. it scatters over several multiples of 2it, while r(x,t) then 
fluctuates around r — 1. The distribution of Fourier components (ftk remains essentially 
unchanged. With dissipative terms and small symmetry breaking H — > the initial con- 
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figuration of fig. 2 slowly relaxes into the line = 0, while for strong symmetry breaking it 
evolves into the quasistable soliton-antisoliton configuration shown in fig. 2 at very late time 
(t = 100). 

The power of some Fourier components |vrfc| 2 (k=2,ll,21,31) of the pion field ir(x,t) (3) 
for the evolution of fig.2 is given in fig. 3 up to t = 15. At time t ~ 10 the final configuration 
shown in fig.2 is essentially established. At earlier times there is no specific enhancement 
for the low-k components. 

To our knowledge the ordering process in 3+1 dimensional 0(4) models with soliton 
stabilizing dynamics has not yet been studied in detail. It might provide interesting features 
of baryon-antibaryon formation through the cooling process of the hot pion gas, as the 
different length scales change as functions of the temperature. But the unexpected and 
anomalous exponents in the scaling behaviour of 1+1-dimensional 0(2) spin systems after 



a quench have recently been connected to topological textures in these models |T6| . 

It may appear as an unpleasant feature that the implementation of a global constraint 
B = n should have such drastic consequences for the time evolution of the local field through 
the requirement that a configuration should never move across the origin. If it is true that 
the gross features of the time evolution are insensitive to the actual value of n, one might 
be tempted to simply ignore this local constraint, i.e. perform the calculation in a cartesian 
parametrization of the full (er-7r) space of the 0(4) model, as it is conventionally done, 
without distinguishing different vacua. This would imply that the total baryon number 
will fluctuate during the time evolution of the field as in a grand canonical ensemble with 
vanishing chemical potential. It is, however, not clear under what conditions that may be a 
legitimate procedure. 

To summarize, if we want to maintain the topological concept of baryon number conser- 
vation in effective chiral theories, we have to observe the angular nature of the chiral field 
and allow for a multitude of physically distinct vacua. The global constraint imposed by the 
actual value of the baryon number may be unimportant. But the fact that the topological 
charge is locally conserved has decisive influence on the dynamics and may invalidade con- 
clusions drawn from the time evolution of chiral fields in the 0(4) model without topological 
constraints. 
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